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ABSTRACT 

We have imaged in CO(2-l) the molecular gas in NGC 1275 (Perseus A), 
the cD galaxy at the center of the Perseus Cluster, at a spatial resolution of 
~1 kpc over a central region of radius ~10 kpc. Per A is known to contain 
~1.3 x 10 10 M Q of molecular gas, which has been proposed to be captured from 
mergers with or ram-pressure stripping of gas-rich galaxies, or accreted from a 
X-ray cooling flow. The molecular gas detected in our image has a total mass of 
~4 x 10 9 M , and for the first time can be seen to be concentrated in three radial 
filaments with lengths ranging from at least 1.1-2.4 kpc all lying in the east-west 
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directions spanning the center of the galaxy to radii of ~8 kpc. The eastern and 
outer western filaments exhibit larger blueshifted velocities with decreasing radii, 
whereas the inner western filament spans the systemic velocity of the galaxy. The 
molecular gas shows no signature of orbital motion, and is therefore unlikely to 
have been captured from gas-rich galaxies. Instead, we are able to reproduce the 
observed kinematics of the two outer filaments as free-fall in the gravitational 
potential of Per A, as would be expected if they originate from a X-ray cooling 
flow. Indeed, all three filaments lie between two prominent X-ray cavities carved 
out by radio jets from Per A, and closely resembles the spatial distribution of 
the coolest X-ray gas in the cluster core. The inferred mass-deposition rate into 
the two outermost filaments alone is roughly 75 M Q yr _1 . This cooling flow can 
provide a nearly continuous supply of molecular gas to fuel the active nucleus in 
Per A. 

Subject headings: galaxies: cooling flow — galaxies: ISM — ISM: molecules - 
radio lines: ISM — galaxies: active — galaxies: individual (Perseus A, 3C 84) 



Introduction 



Rich groups and clusters of galaxies are immersed in hot (^10 7 -10 8 K) X-ray-emitting 
gas that constitutes a large fraction of their baryonic mass (|Sarazinlll986l ). Radiative cooling 
of this gas, by virtue of its X-ray emission, s hould result in an inflow of relatively cool gas t o 
the cluster center; i.e., a X-ray cooling flow (ICowie fc Binneyill977l : iFabian fc Nulsenlll977l ). 
In the absence of any reheating, this flow is inferred to have mass-deposition rates ranging 
10 M yr _1 to ~1000 M yr _1 depending on the given cluster, and over a Hubble 



from 
time 

cD (giant elliptical) galaxy (e.g., reviews by lFabianlll994j : lAllenll2000l ). 



10 10 yr) should deposit ~10 n M t o ~10 13 M m of relatively cool gas in the central 



Numerous searches have been made for the predicted mass sink from X-ray cooling 
flows. These searches cover a broad range of temperatur es, from ionized gas a t temperatures 
of ~10 5,5 K (e.g., through the OVI line in the ultraviolet: iBregman et al.ll2006l. and references 
there in) down to ~10 4 K (e.g., Ha, [Nil] and [SII] lines in the optical; I Crawford &: Fabian 
19921 . and re ferences therein), t o neutral gas at temperatures < 10 3 K in atomic hydrogen 
(21-cm line; lO'Dea et al1ll992l. and references therein), hot molecular hydrogen (near-IR 
vibrational lines; Edge et al.l 120021 . and references therein), as well as cool m olecular hy- 



drogen (traced using ro tational lines of carbon monoxide at mm wavelengths; lEdgd 12001 



Salome fc Combed 120031 . and references therein). None of these searches have found gas at 
quantities anywhere near the abovementioned predicted levels (prior to the era of the XMM- 
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Newton satellite as explained belo w) in putative cooling-flow clusters (e.g., review by [Fabian 
19941 : iMathews fc Brighentill2003l ). 



Over the past few years, observations of hot gas in galaxy clusters with the XMM- 
Newton and Chandra X-ray observatories have radically altered our concept of X-ray cool- 
ing flows. Spectra taken with XMM-Newton at high energy resolution and sensitivity in- 
deed show a decrease in gas temperature towards the cluster center, but fail to reveal 
any detectable gas at temperatures below abo ut one-third the bulk ambient temperature 
( jPeterson et al.l 120031 : iPeterson fc Fabian! 120061 ) . These results provide no evidence for a 
cooling flow, and imply that any gas cooling to temperatures < 1-3 x 10 7 K has a mass- 
deposition rate at least ten times lower than that previously inferred. Furthermore, images 
taken with Chandra at high angular resolutions reveal that the X-ray gas at the centers of 



putative cooling-flow clusters is usually disturbed (e.g.. iFabian et al.ll200CH : McNamara et al. 



2000: Blanton et al. 2001: McNamara et al. 2001: Mazzotta et al 



2002). X-ray cavities are 



commonly seen, spatially coinc ident with radio jets from the central galaxy where detectable 
(IDunn. Fabian, fc Taylorll2005l . and references therein), indicating that energetic particles in- 
jected from an active galactic nucleus (AGN) are responsible for the observed disturbances. 
The work needed to inflate these X-ray cavities implies that the energy output of AGNs is 
more than sufficient to balance radiative losses from the X-ray gas, and can theref ore greatly 
diminish if not quench the cooling flow (jChurazov et al.l 120021 : iBirzan et al.l 120041 ) . 



Against this backdrop, in the last few years searches targeting putative strong cooling- 
flow clusters have become relatively successful at detecting massive amounts of co ol molecular 



hydrogen gas as traced in carbon monoxide (CO) at the centers of these clusters ( lEdgel 12001 



Salome fc Combed I2003T). About twenty such examples are now known, whereas earlier 



searches fe.g.. iGrabelskv fc Ulmer 



19941 : iBraine fc Dupraz 



1994 



cluster, the Perseus Cluster ( iLazareff et al. 



1990 



O'Dea et al. 



McNamara fc Jaffd Il99dl : lAntonucci &: Barvainis 
994) detected CO gas at the center of only one 



19891 : iMirabel et al.l Il989l ) . The earlier searches 



targeted relatively nearby clusters which, in the light of the downward revisions in mass- 
deposition rates, are now known to mostly have very weak if any coolin g flows. On the other 
hand, the more recent surveys (jEdge 2001 : Salome &: Combes 2003 ) target more distant 
clusters with higher mass-deposition rates as inferred in the traditional picture of X-ray 
cooling flows. Follow-up imaging where performed reveals that the molecula r gas is usually 
locat e d within a region of < 20 kpc coincident with the central galaxy (jEdge fc Frayer 



20031 : ISalome fc Combed 120041 ) . The inferred mass of molecular hydrogen gas spans the 



range ~10 9 -10 n M , still much lower than that predicted in the traditional picture of X-ray 
cooling flows, but closer now to the revised picture given the current upper limits in mass- 
deposition rates. The amount of cool molecular gas far outweighs all other gas components 
detected at temperatures below ~10 6 K combined. These results leave open the possibility 
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that X-ray cooling flows may deposit what is still a massive amount of cool gas compared 
with that normally present in galaxies (especially elliptical or cD galaxies). 

A cooling flow, however, is not the only possible source of molecular or other relatively 
cool gas in cD galaxies. Ram-pressure stripping of the interstellar medium of gas-rich galax- 
ies by the hot intracluster medium, or cannibalisms of gas-rich gala xies by the cD galaxy, 
provide plausible if n ot more attractive alternatives (e.g., reviews by Mathews fc Brighenti 
20031 ; iBregmanl 120041 ) . Unlike cooling flo ws, ram-pressure stripp ing is firmly established to 
occur in cluster galaxies (e.g., review bv lBoselli fc GavazzilboOffl. and l ikewise mergers and 
cannibalisms in field galaxies (e.g., review by lBarnes fc Hernquistlll992l ). Indeed, cD galax- 
ies presumably attained their present-day mass through repeated cannibalisms of cluster 
galaxies, or repeated merging of subclusters and galaxy groups before or during cluster viri- 
alization (e.g., see discussion in IPimbblet. Roseboom. fc Doyld 120061 ). In support of these 
arguments, a number of elliptical galaxies in the field, members of groups not recognized 
to have X-ray cooling flows, or those residing in but not at the centers of clusters have 
bee n detected in CO with inferred molecular hydro gen gas masses of up to several ~10 9 M Q 



Knapp fc Rupenl Il996l : iLim et al.l l2000l . I2004J . and references therein) . At least one of 



(see 

these galaxies, NGC 759, resides in but not at the center of a putative cooling-flow cluster, 
Abell 262. NGC 759 has an inferred mass in molecular hydrogen gas of 2.4 x 10 9 M Q , at 
least two orders of magnitude higher than the uppe r limit placed on the m ass of molecular 
hydrogen gas in the central cD galaxy, NGC 708 (IBraine fc Duprazl 119941 ) . These results 
clearly demonstrate that a cooling flow is not the only mechanism capable of depositing 
large quantities of molecular gas in cluster elliptical galaxies. 

At the present time, the available observations do not allow us to confidently discrim- 
inate between the different possible origins for the molecular (or other relatively cool) gas 
in cD galaxies. Yet, naively, such a distinction should be possible simply from the observed 
spatial distribution and kinematics of the molecular gas. In the absence of any disturbances, 
the molecular gas deposited by a X-ray cooling flow should simply flow radially inwards along 
the gravitational potential of the cD galaxy. On the other hand, gas captured from a gas- 
rich galaxy is expected to orbit (and hence exhibit a global rotation) about the cD galaxy, 
reflecting the orbit (or passage) of the gas-rich galaxy around (or through) the cluster center 
or cD galaxy. Here, we present a high angular-resolution CO observation of the cD galaxy 
Perseus A (Per A; NGC 1275) that reveals for the first time an ordered spatial-kinematic 
pattern in its molecular gas. Located at a comparatively close distance of ~74 Mpc (redshift 
z = 0.01756, and assuming H = 70 km s _1 Mpc -1 and Q = 1), Per A lies at the center of 
the Perseus Cluster (Abell 426), the X-ray-brightest cluster in the sky. The radiative cooling 
time in the inner few tens of kpc of the Perseus cluster is less than ~10 9 yr, leading to a 



mass-deposition rate of a few 100 M Q yr 1 in the absence of any heat input (jFabian!ll994 
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Fabian et al.l 12006c iDunn fc Fabian! 120061 ). the highest of any cluster within its distance. 



For over a decade, Perseus A was the only cD galaxy in a putative cooling-flow cluster 
known to contain molecular gas, which together with its other peculiarities and relative 
proximity have made it the subject of inten sive studies. This galaxy has been mapped 
with single-dish telescopes in the J — >■ 1 — (IReuter et al.l Il993l ) , J — >■ 2 — 1 (IReuter et al. 



1993c bridges k Irwi3 Il998c ISalome et all l2006h. and J -> 3-2 (bridges fe Irwinl Il998h 



transitions of CO. The most recent map in CO(2-l) by Salome et al. ( 20061 ). which has 
the highest angular resolution and covers the widest area, shows that the molecular gas 
is elongated approximately east-west across the center of the galaxy. In agreement with 
previous single-dish maps, the molecular gas shows no apparent ordered spatial-kinematic 
pattern, with predominantly blueshifted velocities both east and west of the nucleus, and 
velocities near the systemic velocity immediately around the nucleus. The inferr ed mass of 
molecular hydrogen gas in Per A is ~1.3 x 10 10 M (Fig. 4 of lSalome et al.l 120061 ). far larger 
than that yet detected in field or non-cD cluster elliptical galaxies in the Local Universe (e.g., 
review by lKnapplll999l ) . The only previous inteferometric map of the mo lecular gas i n Per A 
was made in CO(1-0) with the Nobeyama Millimeter Array (NMA) by ilnoue et al.l (119961 ) . 
T his map differs in a number of important respects to single-dish maps (see discussion 



m iSalome et al.l 120061 ). and bears little resemblence to that shown here in CO(2-l). We 
suspect that the CO(1-0) interferometric map may be compromised by inaccurate bandpass 
calibration, as explained in the next section. The only other published attempt to image the 
molecular gas in Per A with an interferometer was made by lBraine et al.l (119951 ) in CO(2-l) 
with the IRAM Plateau de Bure Interferometer (PdBI), where virtually all the emission was 
resolved out. 

Readers interested in the technical details of our observations and data reduction should 
now proceed to §2. Those interested only in the results can skip ahead to §3. In §4 we 
compare the observed properties of the molecular gas with predictions for capture from gas- 
rich galaxies, and in §5 with predictions for accretion from a X-ray cooling flow. In §6, we 
discuss the implications of our results for the cooling flow in the Perseus Cluster. Readers 
interested only in a concise summary of our results and interpretation can proceed directly 
to 87. 



2. Observations and Data Reduction 



Mapping of line emission from Per A at radio wavelengths is made difficult by the much 
stronger continuum emissio n of the galaxy. Ex tended jets and lobes dominate the continuum 
at centimeter wavelengths (jPedlar et al.lll990l ). whereas a compact nuclear source dominates 
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the continuum at millimeter wavelengths (IBraine et al.lll995l ; llnoue et al.lll996l ). In such a 
situation, accurate measurements of the spectral response (passband) of each antenna and 
its variation with time are required. Any measurement errors introduce a frequency and 
time dependence in the point spread function of the telescope, which when convolved with 
the continuum source leave spatially-extended and frequency-dependent artifacts in the map 
that can be mistaken for or alter the line emission. We suspect that such artifacts may be 
responsible for the differences seen between the interferometric map of llnoue et al.l (119961 ) in 
CO(l-O) (rest wavelength at 2.6 mm) and the single-dish maps mentioned in §1, as well as 
the map made here. 

The continuum emission of Per A is weaker at shorter wavelengths, and hence the re- 
quired measurement accuracy and stability of antenna passbands become correspondingly 
less stringent at higher transitions of CO. We have therefore chosen to map Per A in CO(2-l) 
at 1.3 mm, using the SubMillimeter Array (SMA). The observations spanned ~5.5 hrs with 
an on-source time of ~4 hrs on Oct 22, and spanned ~10 hrs with an on-source time of 
~7 hrs on Oct 24. The weather was relatively poor on Oct 22, and so we report here only 
the observation on Oct 24. The SMA was in its compact configuration, and with 6 anten- 
nas operating provided projected baselines over the observation ranging from about 8-68 m 
(6-52 kA). We observed the CO(2-l) line in the upper sideband of the double-sideband re- 
ceivers. The correlator was configured to provide twenty-four independent spectral windows 
("chunks") in each sideband. Each window was divided into 128 channels with a channel 
width of 0.8125 MHz (i.e., ~1.1 km s _1 ), giving a total frequency coverage in each sideband 
(with a small overlap between windows at their edges) of 1.97 GHz (i.e., ~2590 km/s). 

To accurately measure and track variations in the antenna passbands with time, we 
observed Mars at hourly intervals over the duration of our observation. We found that 
the antenna passbands changed only slowly with time, and at a level smaller than that 
required to keep artifacts below the level of random noise fluctuations. We have therefore 
averaged all our measurements of Mars to derive a precise passband solution for each antenna. 
Callisto served as the absolute flux calibrator. We used the strong (~3.4 Jy) and compact 
(unresolved) nuclear continuum source in Per A to derive the complex gains (amplitude and 
phase) of each antenna as a function of time, solving for the phase at 3 min intervals and 
amplitude at 20 min intervals. 

The full-width half-maximum (FWHM) of the SMA primary beam at a frequency of 
226.6 GHz in CO(2-l) is 55", corresponding (with 1" = 360 pc in th e assu med Cosmology) 
to a circle of radius ~10 kpc. The CO(2-l) map of ISalome et al.l (120061 ) shows that the 
molecular gas extends to radii of ~14 kpc both east and west of the nucleus, and hence our 
map spans only the inner central extent of the molecular gas. Our observation was centered 
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on the nucleus of Per A, located at Right Ascension = 03 h 19 m 48!16 and Declination = 
41° 30" 42'.'1. We attained a synthesized beam (with natural weighting of the visibilities) of 
3'/0 x 2'.'7 (1.1 x 1.0 kpc), with major axis at a position angle of 39°. 5. To make maps in CO (2- 
1), we first derived the continuum emission from the line-free channels and subtracted this 
emission from the measured visibilities. We then averaged the data in velocity channels of 
20 km s _1 , and deconvolved the point spread function of the telescope from the image in each 
channel. To make the corresponding maps in total intensity (zeroth moment) and intensity- 
weighted mean velocity (first moment), we first convolved the channel maps with a triangular 
function (i.e., Hanning smoothing) having a full- width at zero intensity of 100 km s _1 , and 
then discarded pixels with intensities less than 1.5a (where a is the root-mean-square noise 
fluctuation in one channel). We then summed the remaining data in the channel maps 
(with no smoothing) to produce the total intensity map, and computed the mean velocity of 
each pixel weighted by its intensity in each channel to produce the intensity-weighted mean 
velocity map. 



Results 



3.1. Spatial distribution and kinematics 



The CO(2-l) channel maps are shown in Figure 1, and the corresponding moment maps 
in Figure 2. In Figure 1, the line intensity in each channel is plotted as contours on a 
grayscale map of the integrated line intensity. In Figure 2, the integrated line intensity 
is plotted as contours on a color-coded map of the intensity-weighted mean velocity. The 
location of the central continuum source, corresponding to the AGN, is indicated by a cross 
in both figures. The velocities are measured with respect t o the systemic heliocentric yelocit y 
of Per A derived in the optical of cz = 5264 ±11 km s _1 (iHuchra. Vogeley. fc Gelled Il999l ). 
where c is the speed of light and z the redshift. For clarity, no primary beam correction 
has been applied to the channel maps shown in Figure 1. On the other hand, to faithfully 
reproduce the overall intensity distribution of the molecular gas, a primary beam correction 
has been made to the integrated line intensity map shown in both Figure 1 and Figure 2. 



As in single-dish maps (IReuter et al.lll993t iBridges & Irwinlll998l ; ISalome et al.l 120061 ). 
the molecular gas in our interferometer map is concentrated primarily in the east-west direc- 
tions. Also as in single dish maps, the emission away from the nucleus on both the east and 
west, as well as north and south, is entirely blueshifted, with emission closest to the nucleus 
redshifted. For the first time, our interferometer map reveals that the molecular gas in the 
central region of radius ~10 kpc is concentrated primarily in three radial filaments, labeled 
in Figure 2 as the inner, western, and eastern filaments in order of increasing projected radii 
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from the nucleus. (As explained below, the inner filament may not actually comprise a single 
integral filament but rather two separate structures.) Although the inner filament appears 
to connect with the western filament, they span very different velocities and have entirely 
different radial velocity gradients. These filaments appear as distinct spatial-kinematic fea- 
tures in the channel maps (Fig. 1), and are clearly separate structures. Similarly, the feature 
labeled El in Figure 2 is discontinuous in velocity from the eastern filament, and is regarded 
here as a separate structure (blob). Another blob labeled E2 can be seen just to the north of 
the eastern filament. An apparently isolated blob labeled Wl lies to the west of the western 
filament, located along a line from center that nearly passes through the inner and western 
filaments. The only features located away from the east-west directions are the isolated 
blobs Nl and SI located to the north and south of center respectively. The spectra of all 
the features labeled in Figure 2 are shown in Figure 3. 

O ur CO(2-l) interfer ometer map bears little resemblance to the CO (1-0) interferometer 



map of llnoue et al.l (119961 ). which as mentioned in §1 also shows significant differences com- 
pared with single-dish maps. First, none of the eastern features seen in Figure 2 are detected 
in the CO (1-0) interferometer map, despite the presence of gas at this location in single-dish 
maps. Second, the inner and western filaments do not appear as distinct features in the 
CO (1-0) map, but coincide with a complex and spatially extended distribution of gas in the 
CO(l-O) map. Third, gas in the CO(l-O) map at the location of the western filament ex- 
hibits both blueshifted and (primarily) redshifted velocities; in our map, the western filament 
exhibits only blueshifted velocities, consistent with that seen at this location in single-dish 
maps. Finally, the CO(l-O) map shows gas towards the south-east of the nucleu s that is not 



seen i n our CO(2-l) map. We conclude that the CO(l-O) interferometer map of llnoue et al. 



(119961 ) is unreliable, and cannot be meaningfully compared with our CO(2-l) map. 



The radial positions and extents, as well as velocities, of all the features labeled in 
Figure 2 are listed in Table 1. To estimate the lengths of the filaments, we measure from 
the channel maps the difference in radial positions of the emission centroids at opposite ends 
in velocity spanned by that filament. Because the emission is spatially resolved in many 
channels, this method provides a conservative estimate of the true length of the filaments, 
but is necessary in particular for the eastern and western filaments as they do not have 
a smoothly-varying intensity profile along their lengths. As measured from the emission 
centroid in the channel maps, the inner filament can be traced from a (projected) radius of 
l'/2 ± 0'.'3 (432 ± 108 pc) where it has a (projected) velocity of 100 km s _1 to a radius of 
4'.' 2 ± 0'/3 (1.5 ± 0.1 kpc) westwards where it has a velocity of —100 km s _1 , thus spanning a 
(projected) length of ~3'.'0 (~1.1 kpc). A two-dimensional Gaussian fit to the inner filament 
in the integrated intensity map of Figure 2 gives a length (at full-width half maximum) 
of 4'.'7 ± 0'.'6 (1.7 ± 0.3 kpc). This filament straddles the systemic velocity of the galaxy, 
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and exhibits a radial velocity gradient with smaller blueshifted velocities or larger redshifted 
velocities with decreasing radii (Fig. 2). 

The western filament can be traced from a radius of 9 / /0 ± 0'/3 (3.2 ±0.1 kpc) where 
it has a velocity of —180 km s _1 to a radius of 15'/7 ± 0'/4 (5.7 ± 0.1 kpc) westwards where 
it has a velocity of —80 km s _1 , thus spanning a length of ~6'/7 (~2.4 kpc). The eastern 
filament can be traced from a radius of 15'.'2 ± 0'.'4 (5.5 ± 0.1 kpc) where it has a velocity 
of —160 km s _1 to a radius of 21'/3 ± 0'/3 (7.7 ± 0.1 kpc) eastwards where it has a velocity 
of —40 km s , thus spanning a length of 6'.'1 (2.2 kpc). Unlike the inner filament, the two 
outer filaments (and all the blobs) exhibit entirely blueshifted velocities, consistent with the 
velocities measured in single-dish maps at their locations. These filaments can be seen to 
exhibit a radial velocity gradient with larger blueshifted velocities at smaller radii, opposite 
to that exhibited by the inner filament. 

In Figure 4, we plot the CO(2-l) intensity of each filament along their radial extents as 
a function of velocity (i.e., position- velocity (PV-) diagrams). For clarity, no primary beam 
correction has been applied to these plots. Both the eastern and western filaments exhibit a 
simple linear change in velocity with radius from the center of the galaxy. The feature at a 
position of 26" and velocity of —80 km s _1 in the PV-diagram of the eastern filament (bottom 
panel of Fig. 4) corresponds to the feature El (Fig. 2), which as mentioned earlier probably 
consitutes a separate structure. The PV-diagram of the inner filament is more complicated: 
the velocity changes linearly with radius at both blueshifted and redshifted velocities, but 
with different velocity gradients on either sides of the systemic velocity. It is possible that the 
redshift and blueshifted sides of the inner filament correspond to separate structures that are 
poorly resolved spatially in our map. Indeed, a close examination of Figure 2 suggests that 
the redshifted and blueshifted portions of this filament lies at somewhat different position 
angles, producing a visible twist in the major axis of this filament between the two portions. 



3.2. Molecular gas mass 



The integrated line intensities, corrected for the primary beam response of the SMA, 
and corresponding molecular hydrogen (H2) gas masses of all the features labeled in Figure 2 
ar e listed in Tab l e 1. T o facilitate a direct comparison with the molecular gas masses derived 



111 



Salome et al.l (120061 ). we adopt like them a brightness temperature in CO (2-1) of 0.7 times 



that in CO(1-0), and the standard C O (1-0) to H2 conversi on ratio (applicable to gas with 



solar metallicity) of 2.3 x 10 cm flSolomon et al.l 119971 ). Altogether, the features have 



a total gas mass of (4.2 ± 0.2) x 10 9 M , about half that detected i n observations wit h 
single-dish telescopes over the same radius of ~1 x 10 10 M (Fig. 4 of ISalome et al.l 120061 ). 
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The three filaments alone have a total gas mass of (3.3 ± 0.2) x 10 9 M , about 80% of that 
contained in all the features. Our interferometric observation is preferentially sensitive to 
relatively bright and compact features, and hence traces regions where the molecular gas is 
most strongly concentrated (i.e., highest filling factor). 

To determine if there is emission external to the features detected in Figure 2 that lies 
below our brightness sensitivity, we have integrated the channel maps over the entire SMA 
primary beam from a velocity of —230 km s _1 to 170 km s" 1 . The integrated line intensity 
determined in this manner is 409 ± 28 Jy km s _1 , corresponding to a molecular gas mass 
of (11.8 ± 0.8) x 10 9 M . This is higher than the total molecular gas mass of the features 
detected in Figure 2, and comparable with the abovementioned molecular gas mass detected 
with single-dish telescopes over the same region. 



3.3. Correspondence with Ha gas 

Per A exhibits a complex H a nebula that extends beyond the optically visible confines 



of the galaxy to nearly 70 kpc (jConselice et al.ll200ll ). The origin and excitation of this 



nebula is very poorly understood. In Figure 5, we plot the integrated CO(2-l) intensity on a 
color map of the inner region of this Ha nebula. All the molecular filaments and blobs lying 
in the east-west directions, but not the blobs lying north and south of center, appear to be 
associated with the Ha nebula. The inner filament closely coincides with a bright western 
extension in Ha from a compact nuclear component. The inner radius of the western filament 
closely coincides with an apparent extension of the bright central Ha feature further west. 
The eastern filament coincides with local brightenings in Ha, as do the nearby blobs El and 
E2. 



Ferruit et al.l (119971 ) have measured the velocity of the Ha-emitting gas within a radius 
of ~10" of the center of Per A, approching to within 2 / . / 6 of its nucleus. Along the inner 
filament, the Ha emission has a velocity of about —25 km s _1 at a radius of 2'.'6 and about 
— 105 km s _1 at a radius of 3'.'5. This radial range coincides with the outer portion of the inner 
filament where the CO (2-1) emission is primarily blueshifted: the Ha velocity compares well 
with the range of CO (2-1) velocities seen here, including the trend of increasingly blueshifted 
velocities towards larger radii. Along the eastern filament, the Ha emission has a velocity of 
about —165 km s _1 at a radius of 5'/7, which coincides quite closely to the inner portion of the 
eastern filament. Here, the CO(2-l) velocity reaches about —180 km s , again comparable 
with that seen in Ha. 
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Conselice et al.l (120011 ) have measured the velocity of the Ha-emitting gas out to a radius 



of nearly 10" on the eastern side and nearly 20" on the western side, but with poorer spatial 
sampling. At and close to the nucleus, the Ha emission is redshifted, like in in our CO(2-l) 
map. Further west at the location of the western filament, the gas is blueshifted, again as 
is seen in our CO(2-l) map. We therefore conclude that the CO(2-l) features are not just 
spatially coincident with bright Ha features but, where the velocity of the latter has been 
measured, also have comparable velocities. 



3.4. Correspondence with X-ray features 



In Figure 6, we overlay (black) contours of the inte grated CO (2- 1 ) inten sity on a X-ray 
image of the center of the Perseus cluster taken from iFabian et al.l (120061 ). We also plot 
in (red) contours the radio e mission of Pe r A at 90 cm, pr oduced by relativis tic jets from 
the central AGN, taken from iPedlar et al.l (119901 ) (see also IFabian et al.l 120021 ). The X-ray 
image shows cavities to the north and south of the nucleus that coincide with the radio jets. 
These cavities are therefore believed to be giant bubbles blown by the radio jets into the 
surrounding X-ray gas. The bulk of the molecular gas — specifically, all the features located 
east and west of center — lies (in projection) between the X-ray cavities. The isolated blob 
Nl lies close (in projection) to a local X-ray brightening projected against the northern X-ray 
cavity, and the blob SI near the inner rim of the southern X-ray cavity. 

In Figure 7, we overlay contours of the integrated CO (2-1) intensity on a map of the 
lowest temperature (0.5 keV, or ~5 x 10 6 K) X-ray ga s dete cted by Chandra towards the 
center of the Perseus cluster, taken from IFabian et al.l (120061 ). The bulk of this relatively 
cool X-ray gas lies in an east- west band (curving north) across the c e nter o f the galaxy, 
located between the two X-ray cavities. As remarked by IFabian et al.l (120061 ). the bulk of 
this relatively cool X-ray gas clos ely resembles the spatial distribution of the Ha gas (see 
also Fig. 12 of IFabian et al.l 120061 ) . although the mass of the X-ray gas at this temperature 
alone is an order of magnitude larger than that of the Ha gas. The bulk of the molecular 
gas detected in our image also has a spatial distribution that closely matches that of this 
relatively cool X-ray gas. 



4. Gas captured from a gas-rich galaxy? 



The origin of the molecular gas in Per A has long been a subject of debate and con- 
troversy. Arguments have been made for mergers with or cannibalisms of gas-rich galaxies, 
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ram-pressure stripping of gas-rich galaxies passing close to or through Per A, and a X-ray 
cooling flow. Here, we examine whether the observed properties of the molecular gas are 
consistent with capture from gas-rich galaxies. 



4.1. Interactions with a foreground infalling galaxy 



Suggestions that the molecular gas in Per A is captured from either ram-pressure strip- 
ping or mergers/cannibalisms of gas-rich galaxies are bolstered by evidence that Per A may 
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the continuum emission of Per A at a velocity around 8200 km s -1 . Despite having a ve- 
locity that is ~3000 km s^ 1 higher than that of Per A, and hence commonly referred to as 
the high- velocity system (HVS), the detecti on of line absorption in both HI and Lyman a, 
as we ll as continuum absorption in X-rays (IFabian et al.ll2000l : iGillmon. Sanders, fc Fabian 
2004 ). implies that the HVS lies in front of Per A. 



Evidence that this galaxy is not just a chance foregound galaxy but in the close vicinity 
of Per A comes from the detection of Ha emission at velocities intermediate between those 



of Per A and the HVS (IFerruit et al.lll997l ). This intermediate- velocity gas may comprise gas 
from the HVS entrained by the hot intracluster X-ray ga s. Furthermore, the ionized gas i n 



the HVS appears to exhibit a complex velocity field fe.g.. lBorosonlll990l: lUnger et al.lll990l ). 
which together with the active star formation in this galaxy (ICaulet et al.lll992l ) suggests 
that the HVS has been disturbed by the intracluster medium or interactions with Per A. 

It is highly unlikely, however, that the molecular gas in Per A was captured from the 
HVS. First, the velocity of the molecular gas is ~3000 km s" 1 lower than the velocity of the 
HVS, and is predominantly blueshifted with respect to the systemic velocity of Per A rather 
than redshifted like the HVS. Second, the HVS has not been detected i n HI emission nor CO, 
placi ng a 3a upper mass limit of 2 x 10 9 M in at omic hydro gen gas (Ivan Gorkom fc Ekers 
19831 ) and 5 x 10 9 M in molecular hydrogen gas (I Jaffd 119871 ). The latter is about a factor 
of four lower than the mass of molecular h ydrogen gas in Per A. Consistent with these argu- 
ments, IGillmon. Sanders, fc Fabian! (12004 ) place a lower limit of ~60 kpc on the separation 
between the HVS and Perseus A, and argue that any collision between the two galaxies still 
lies in the future. 
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4.2. Previous encounters with gas-rich galaxies 



If captured, the molecular gas in Per A more likely originated from previous encounters 
with other gas- rich galaxies. Gas captured in this manner should exhibit a transverse (i.e., 
perpendicular to radial) as well as any radial spatial-kinematic components following the 
path of the interacting galaxy through or decaying orbit of the cannibalized galaxy in the 
cD galaxy. In both cases, one might expect the captured gas to orbit and hence exhibit a 
global rotation about the cD galaxy. 

Where imaged and well resolved, the CO gas in elliptical galaxies has until now been 



invariably found to be distributed in rotating disks (e.g., see lYoungl 120051 . and references 
therein), suggesting that their gas was indeed captured from gas-rich galaxies. Where only 
CO spectra (at good signal-to-noise ra tio) are available, many exhibit double-horned line 
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profiles indicative of rotating disks (e.g. , iMazzarella et al.lll993l ; IWiklind. Combes, fc Henkel 



kinematic components nor any ordered global motion indicative of rotation. Instead, the 
most prominent features that we detected in Per A — the three filaments described in §3 - 
exhibit only radial spatial-kinematic components. If each filament is attributed to capture 
from an individual gas-rich galaxy passing through Per A, this would require three separate 
events with each gas-rich galaxy taking a projected path that crosses the center of Per A. 
Although this could be the future fate of the HVS, there is no independent evidence that 
three such highly unlikely events have occurred in the relatively recent past. 



5. Gas accreted from a X-ray cooling flow? 



The excellent spatial correspondence between the molecu lar gas and the rela tively cool 
X-ray gas (Fig. 7), but not the warmer X-ray gas (Fig. 12 of iFabian et al.l 120061 ). in Per A 
provides circumstantial evidence that the molecular gas originates from a cooling flow. In the 
simplest possible situation imaginable, gas accreted from such a cooling flow should simply 
flow radially inwards along the gravitational potential of the cD galaxy. The gas should 
therefore exhibit only radial spatial-kinematic components, as is indeed observed for the 
three filaments. Such a signature alone, however, is not a unique indicator of radial inflows. 
In principle, gas undergoing radial outflows also can exhibit the same spatial-kinematic 
pattern. We now examine in greater detail the properties of these filaments to see if we can 
distinguish between inflows or outflows. 
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5.1. Correspondence with dust 



5.1.1. Outer molecular filaments and blobs 

The eastern and western filaments, as well as the other blobs of molecular gas, exhibit 
only blueshifted velocities. This implies inflow if the filaments and blobs are located on the 
far (back) side of Per A, or outflows if located on the near (front) side. To deduce which is 
the case, we have compared the spatial distribu t ions of these filaments and blobs to optica l 



images showing dust (e.g., ICarlson et al.l Il998l ; iKeel &: White! l200ll ; IConselice et al.ll200ll ) 



seen in silhouette against starlight from Per A. If these filaments and blobs are located on 
the near side of the galaxy, they ought to show a good correspondence with the silhouette 
dust. On the contrary, if these filaments and blobs are located on the far side of the galaxy 
and at a large inclination to the plane of the sky, they may not show up strongly if at all in 
silhouette dust. Examples of s ilhoue t te dust seen in elliptical gala x ies with the HST-W FPC2 
can be found inlMartel et all Jl999h . IVerdoes Kleiin etaD Jl999h . Ide Koff et aD (boooh . and 
Tran et al.l (120011 ). Inferred dust masses as small as 10-100 M are detectable, which for a 
ratio in mass of molecular hydrogen gas to dust of about 100 (as in our Galaxy) corresponds 
to molecular gas masses as small as just 10 3 -10 M . 

In Figure 8, we overlay the integr ated CO(2-l) intensity on the HST-WFPC2 image of 
Per A taken from lCarlson et al.l (119981 ). An especially prominent and complicated network of 
dust can be seen in silhouette towards the north of the nucleus, extending both east and even 
furth er west. Much if not all of this dust likely belongs to the foreground infalling galaxy 
(e.g., IKeel fc White! l200ll ; IConselice et al.ll200ll ). First, the overall spatial distribution of 
this dust network closely resembles that of the Ha gas in the HVS, but not that in Per A 
itself (e.g., compare th e silhouette dust in Fig. 7 with the individual Ha images of the HVS 
and Per A in Fig. 3 of ICaulet et al.lll992l ). In addition, the morphology of this dust closely 
matches that of gas detected in X -ray absorption against Per A, attributed to the HVS 
( iGillmon. Sanders, fc Fabian! 120041 ) . 



Whether belonging to the HVS or Per A, or portions in both, neither the eastern nor 
western filaments, nor the blobs El, E2, Wl, and Nl coincide with any of the silhouette 
dust (SI is located beyond the region spanned by the HST-WFPC2 image). Specifically, 
the complicated network of dust does not extend to the location of the eastern filament and 
blobs, where the light distribution appears to be smooth. Dust in silhouette is visible to the 
north, but not at the location, of the western filament and Wl. The light distribution at and 
immediately around the western filament and Wl appears to be quite smooth, although its 
proximity to the prominent silhouette dust does make any weaker silhouette features difficult 
to discern. The lack of correspondence with silhouette dust suggests that both the eastern 
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and western filaments, as well as the other blobs of molecular gas, are located on the far side 
of the galaxy and at a relatively large inclination to the plane of the sky. Their blueshifted 
velocities therefore imply inflow. We disregard the possibility that these filaments and blobs 
lack dust because they formed from me tal-poor gas; on the contrary, the iron abu ndance 
of the X-ray gas is at least ~0.5 solar (ISchmidt et al.l 120021 ; IChurazov et al.l 120021 ) within 
~100 kpc. 



5.1.2. Inner molecular filament 

To permit a more detailed comparison between the inner molecular filament and the 
silhouette dust, in Figure 9 we overlay the integrated CO(2-l) intensity o n just the inner 



portion of HST-WFPC2 image of Per A taken from Keel & White! (1200 if ). The latter is 



a transmission map, representing the residual intensities after subtracting smooth model 



profil es for both the foreground and background galaxy (for more details, see iKeel fc White 



20011 ). Brighter regions therefore indicate intensities above the smooth background, and 



darker regions intensities below the smooth background (i.e., dust extinction). 

Confirming what was seen in Figure 8, neither the western molecular filament nor blob 
Wl obviously coincide with the silhouette dust. There is a patch of dust at the innermost 
radii of the western filament where the CO(2-l) emission is strongest, but this dust patch 
does not trace the length of the filament and hence is likely coincidental. There is patchy 
filamentary dust at and around the location of the inner molecular filament, although once 
again this dust does not obviously trace out the filament. If the dust seen here also belongs 
to the HVS, this would make any silhouette dust close to the nucleus of Per A difficult 
to discern. Thus, at the present time, we cannot be certain whether the inner filament is 
associated with any silhouette dust, and hence whether this filament is located on the near 
or far side of Per A. 



5.2. Correspondence with radio jets 

The relative locations of the molecular features and radio jets from Per A were shown 
earlier in Figure 6. If the molecular filaments we observe correspond to gas entrained by the 
radio jets, they should lie in projection against, or at the edges of, these jets. In addition, the 
filaments should have their long axes aligned approximately north-south, along the axis of 
the radio jets. Instead, the molecular filaments have their long axes aligned approximately 
east-west, closely orthogonal to the axis of the radio jets. Furthermore, all the filaments 
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lie neatly between the X-ray cavities, in regions where there is presumably no energetic 
particles to excavate the X-ray gas. Thus, the molecular filaments are unlikely to comprise 
gas entrained by the radio jets or other energetic particles. The same argument applies to 
the blobs El, E2, and Wl. 

On the other hand, the blobs Nl and SI do lie in projection against the radio jets. SI 
is spatially resolved, and appears to be elongated in an approximately west-west direction, 
closely parallel with the western filaments. If entrained by the radio jet, one would expect 
this blob to be elongated approximately north-south, along the axis of the radio jet. 



5.3. Free-fall in gravitational potential of Per A 

Both the eastern and western filaments exhibit a linear increase in blueshifted velocities 
with decreasing radii. Here, we examine whether such a motion can be produced by gas in 
free-fall along the gravitational potential of Per A. 



5.3.1. Model for gravitational potential 
To m o del th e gravitational potential of Per A, we use an analytical model proposed by 



Hernquistl (119901 ) that has been found to reproduce the projected surface brightness profile 
of elliptical galaxies parameterised only by their mass and effective radius. In this scheme, 
the mass density function of radius is given by 

Ma 1 , , 

p[r) 



2-7T r (r + a) A 

where M is the total mass of the galaxy, r the radius, and a a scale length. The scale length 
a is related to the half-mass radius, ri/2, which contains half of the mass of the galaxy by 

n /2 = (1 + V2)a . (2) 
The gravitational potential is then given by 

#M = -— , (3) 



where G is the gravitational constant. The velocity, v(r), of an object in free-fall can be 
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computed from the relation 

v(r) 2 GM v{r f GM 



2 (r + a) 2 (r + a) 
where r c is the initial radius at which the object is released and v(r Q ) its initial velocity. 



(4) 



Taking values from ISmith. Heckman. &: Illingworthl (119901 ) but scaled to the Cosmol- 
ogy used here, Per A has an optical V-band luminosity of 1.1 x 10 11 L Q and an esti- 
mated mass-to-light ratio in this band of M/L ~ 7.4. The inferred mass of the galaxy 
is therefore 8.3 x 10 11 M ; given the uncertainty in M/L (which, moreover, may change 
with radius), there is probably cons iderable uncertainty in the estimated galaxy mass. 



Smith. Heckman. fc Illingworthl (Il990l ) quote an effective radius, R e , for Per A of 12.4 kpc. 



The effective radius is related to the scale length a by 

R e « 1.8153a , (5) 

and so a ~ 6.8 kpc. 



5.3.2. Predicted velocity profile 
Outer filaments 

The X-ray gas is very nearly in hydrostatic equilibrium, and so a given parcel of cooling 
gas is initialy at rest with respect to the cluster. In Figure 10, we plot from Eq. (4) how 
an object, which is initially at rest, changes in velocity with radius when dropped from two 
different initial radii, at 8.5 kpc (upper panel) and 12.2 kpc (lower panel). (The choice of 
initial radii is explained in the next paragraph.) For comparative purposes, but also because 
of the considerable uncertainty in the estimated galaxy mass, we show the results for two 
different masses for Per A of 3.4 x 10 11 M (solid curve) and the value estimated above of 
8.3 x 10 11 M (dashed curve). As can be seen, beyond a distance of ~1 kpc from which 
the object is dropped, the velocity increases essentially linearly with decreasing radius, as 
is observed for both the eastern and western filaments, irrespective of the assumed galaxy 
mass. The challenge here is to reproduce not just the observed velocity gradient, but also 
the actual velocities spanned by the two filaments. 



With a suitable choice for the inclination of the filaments as the only free parameter, we 
find that this simple model can reproduce the measured range of velocities of both the eastern 
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and western filaments provided we adopt the lower value for the galaxy mass of 3.4 x 10 11 M . 
In Figure 10, we overplot the PV-diagrams of the western and eastern filaments shown in 
Figure 4 but deprojected in both velocity and radius for inclinations that give the best 
model fit. As discussed in §5.1.1, these two filaments must be located at a relatively large 
inclination to the plane of the sky (and on the far side of Per A) for dust in silhouette not to 
be seen associated with these filaments. We find relatively good fits to the predicted curves 
for an inclination of about 40° for the western filament and 47° for the eastern filament. 
Assuming larger inclinations would result in a larger velocity gradient consistent with a 
higher galaxy mass (dashed curves in Fig. 9), but lower deprojected velocities placing these 
filaments on the part of the curve where the velocity increases non-linearly with decreasing 
radius. Assuming smaller inclinations, on the other hand, would result in a smaller velocity 
gradient inconsistent with the model (unless one is prepared to accept an even lower galaxy 
mass). The choice of initial radii mentioned in the previous paragraph is therefore determined 
by our model once a galaxy mass is adopted and the best fit inclinations for the filaments 
determined. 

For an inclination of 40°, the western filament (upper panel) spans deprojected radii of 
4.2-7.4 kpc. In this model, the outermost tip of the filament lies just 1.1 kpc away from 
the location where it first began its free fall in the gravitational potential of Per A. The 
deprojected velocity of this filament at its outermost radius is about —125 km s _1 compared 
with the computed model value of about —121 km s _1 , and at its innermost radius about 
—280 km s" 1 compared with the computed model value of about —273 km s _1 . For an 
inclination of 47°, the eastern filament (lower panel) spans deprojected radii of 8.1-11.3 kpc, 
with its outermost tip lying just 0.9 kpc away from where it began its free fall (i.e., at a 
radius of 12.2 kpc). The deprojected velocity of this filament at its outermost radius is 
about —55 km s _1 compared with the computed model value of about —90 km s -1 , and at 
its innermost radius about —190 km s _1 compared with the computed model value of about 
—206 km s _1 . We regard the agreement between this simple model and our observations 
to be satisfactory given that we have ignored a number of effects, primarily the outward 
pressure the X-ray gas exerts on the infalling molecular gas. The latter should cause an 
infalling object to increase more slowly in velocity with decreasing radius (i.e., resulting in 
a shallower velocity gradient) than predicted by our simple model, thus requiring a larger 
galaxy mass to produce the observed velocity gradient than assumed in our best-fit model. 
This would bring our results into better agreement with model predictions for the inferred 
larger galaxy mass for Per A of 8.3 x 10 11 M Q . 



The scenario we have in mind then is the following. The CO(2-l) transition is excited 
at a critical density in molecular hydrogen gas of ~10 3 cm -3 . Fro m an analysis of the line 
ratios in CO(2-l) and CO(1-0) based on single-dish observations, iBridges fc Irwin! (119981 ) 
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also infer a density for the molecular hydrogen gas of ~10 3 cm -3 , and a temperature of 
20-30 K. By comparison, the density of th e X-ray gas at radii o f 5-10 kpc from the center 
of Per A lies in the range 0.05-0.1 cm -3 (jChurazov et all 120021 ; iFabian et al.ll2006l ). The 
western filament corresponds to molecular gas that has condensed from X-ray gas located 
at a distance of ~8.5 kpc from the center of Per A. The cooling gas decouples from the 
X-ray gas, and then free falls in the gravitational potential of Per A. By the time it has 
travelled 1.1 kpc from its initial location, taking about 18.2 Myr, the cooling infalling gas 
has already been (largely) converted to molecular gas and has attained a density sufficiently 
high and temperature sufficiently low to be detectable in CO(2-l). Because the filament is 
relatively dense, it is unlikely to be (significantly) disrupted during its passage through the 
more diffuse X-ray gas. Furthermore, its linear appearence suggests that the filament is not 
significantly perturbed by turbulence or any other disturbances. The same picture applies 
for the eastern filament, except that in this case the molecular gas condensed from X-ray 
gas initially located at a distance of ~12.2 kpc from the center of Per A. The time it takes 
for the cooling infalling gas to travel 0.9 kpc from its initial location to the outermost tip of 
the eastern filament is about 21.1 Myr. 



Inner filament 

The measured properties of the inner filament, on the other hand, are difficult to explain 
with a simple model. Consider the situation where this filament also corresponds to radially- 
infalling molecular gas, as we have shown is likely the case for the two outer filaments. The 
blueshifted portion of the filament would then have to lie on the far side of the galaxy, 
and the redshifted portion on the near side. As mentioned in §3, these two portions of the 
filament may actually be separate features that are poorly resolved spatially in our map. 
Unlike the two outer filaments, however, the blueshifted portion of the inner filament exhibit 
smaller blueshifted velocities at smaller radii (see top panel of Fig. 4). It may be that, close 
to the center of the galaxy, the infalling molecular gas is being retarded by the high pressure 
of the X-ray gas, such that the blueshifted portion of the filament experiences an inward 
decrease in velocity. The redshifted portion of the filament, on the other hand, exhibits 
larger redshifted velocities at smaller radii, implying an inward increase in velocity. 

Alternatively, and perhaps more plausibly, this filament may correspond to infalling 
molecular gas that is settling into the gravitational potential well of Per A. Specifically, the 



gas may be set t 



ing i nto orbits to form a dynamically unrelaxed disk around the nucleus. 



Heckman et all (119851 ) tentatively infer an approximately constant stellar orbital velocity of 
~50 km s _1 over radii of ~5"-25" in Per A. It is reasonable to expect that as one approaches 
the center of Per A, the X-ray gas increasingly partakes in the overall galaxy rotation. In 
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such a case the molecular gas deposited by the X-ray cooling flow should also possess angular 
momentum, and cannot directly flow into the center but must eventually form a rotationally- 
supported structure (i.e., disk) around the nucleus. 

In support of this picture, hot (~1000 K) molecular hydrogen gas has been detected in 

the n ear-IR vibrational line of H2 v = 1—0 S(l) around the nucleus of Per A (IWilman. Edge, fc Johnstone 



20051 ). This gas appears to be distributed in a disk-like structure oriented east- west with a ra- 



dius of ~50 pc, and therefore shares the same morphological orientation as the inner filament. 
The eastern side of the H 2 v = 1 — S(l) disk is redshifted in velocity by ~100 km s -1 and 
the western side blueshifted by ~100 km s -1 , therefore spanning the same range of velocities 
and also having the same sense in velocity gradient as the inner filament. The inner filament 
may therefore trace the outer, cooler, dynamically unrelaxed portion of the inner, warmer, 
and more dynamically relaxed molecular hydrogen disk, and is responsible for despositing 
gas into this circumnuclear molecular disk. 



6. Implications for Cooling Flow 
6.1. Mass-deposition rate 

6.1.1. Outer molecular filaments 

Given that the eastern and western filaments exhibit a linear dependence in velocity 
with radius, we can compute their dynamical ages from the relationship 

m Tmax ^min , / Umax \ / n \ 

Tdyn = —. r In I , (6) 

where v{r max ) is the deprojected velocity at the maximum deprojected radial extent r max , and 
v(r min ) is the deprojected velocity at the minimum deprojected radial extent r min . For the 
values inferred in §5.3.2, we derive a dynamical age of ~16 Myr for the western filament and 
~29 Myr for the eastern filament. The analytical model described in Sect 5.3.2 (where the 
predicted velocity changes close to but not perfectly linearly with radius) gives a dynamical 
age of 21.7 Myr for the western filament and 16.5 Myr for the eastern filament. Changing the 
inferred inclinations of the filaments by ±20° only changes their dynamical age by a factor 
of ±2 (the dynamical age is proportional to tan(i), where % is the inclination of the filaments 
to the plane of the sky). 

Given that the western filament has a molecular gas mass of ~8.5 x 10 8 M and the 
eastern filament ~6.1 x 10 s M as listed in Table 1, the total mass-deposition rate into 
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these two filaments is therefore ~75 M yr -1 . This value should be considered as a very 
rough estimate that is probably accurate to only a factor of a few. It would be instructive 
to compare the inferred mass-deposition rate into these two molecular filaments with the 
upper limit in mass-deposition rate inferred for the X-ray gas at the same radius. At the 
present time, however, no such information for the Perseus Cluster has been published. 
Instead, below we compare the inferred mass-deposition rate with that inferred for other gas 
components below X-ray temperatures. 



6.1.2. Other gas components 



Bregman et al.l (120061 ) have detected O VI lines tracing ionized gas at temperatures of 
~10 5,5 K in a 30" square region centered on Per A. They attribute these lines to gas from 
a X-ray cooling flow, and infer a mass-deposition rate of 32 ± 5 M Q yr _1 within a radius 
of ~6 kpc. The inferred mass-deposition rate is roughly comparable with that inferred for 
the band of X-ray gas at ~5 x 10 6 K located betwee n the X-ray cavities, provided this 
gas originates from a cooling flow (IFabian et al.l 120061 ). If responsible in its entirety for 
producing (at the present mass- deposit ion rate) the inner molecular filament, which has a 
mass of 1.8 x 10 9 M (Table 1), the cooling flow in the inner ~6 kpc regio n of the cluster 
must have lasted at least ~56 Myr. The mass- deposit ion rate inferred by iBregman et al. 



( 120061 ) is, coincidentally, roughly comparable with the mass-deposition rate we infer for the 
two outer molecular filaments. These filaments, however, lie beyond the region spanned by 
the observed O VI lines. In addition, there is molecular gas lying even further out than the 
region we have mapped (and which also may originate from a cooling flow), suggesting that 
the X-ray cooling flow in the Perseus cluster may have a significantly higher mass-deposition 
rate than that inferred from the O VI lines alone. 

The close spatial-kinematic correspondence between the cool molecular gas detected 
here and ionized gas at ~10 4 K traced in Ha suggests that at least a portion of the latter 
could also originate from a X-ray cooling flow. The mass in Ha gas over its entire visible 
extent is only ~10 M Q ( jHeckman et al.lll989l ). two orders of magnitude below the mass 
in molecular gas of the filaments reported here, and three orders of magnitude below the 
total mass of molecular gas in Per A. The bulk of the gas from the X-ray cooling flow must 
therefore spend relatively little time at temperatures of ~10 4 K. 

One of the outstanding challenges is to understand the time evolution of the cooling flow 
as it transforms from X-ray-emitting plama into molecular gas. At a radius of ~10 kpc where 
we infer the two outer molecular filaments to first cool fr om the X-ray gas, the azimuthally- 
averaged cooling time of the X-ray gas is ~500 Myr (IDunn & Fabian! 120061 ). Although 
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the cooling flow must spend relatively little time at temperatures of ~10 4 K, it may well 
spend much longer times at other temperatures before condensing into molecular hydrogen 
gas. Specifically, the time it spends at temperatures above those traced in CO(2-l) may 
be considerably longer than the inferred free-fall time of just ~20 Myr to reach the outer 
tips of the molecular filaments. If the cooling gas remains coupled to the X-ray gas by 
magnetic fields, an d strong magnetic fields have been inferred in the central region of Per A 
( Taylor et alil2006l ). then this gas may not start to free-fall until atomic hydrogen has (almost 
fully) recombined at temperature below ~4000 K. 



6.2. AGN fueling and feedback 



Molecular gas accreted from a X-ray cooling flow is an excellent candidate for fueling 
the AGN observed in many cD galaxies lying in rich clusters. In the case of Per A, the 
circumnuclear hot molecular hydrogen disk detected by IWilman. Edge, fc Johnstone! (120051 ) 
may be responsible for fueling the AGN; this disk in turn may have been accreted from or is 
replenished by the inner cool molecular filament observed here (§5.3.2). The powerful radio 
jets produced by the AGN rehea ts the surrounding X - ray gas, thereby di minishing; if not 



quenching the cooling flow (e.g., IChurazov et al.l 120021 ; iBirzan et al.ll2004l ). If the cooling 



flow is severely disrupted, the AGN may eventually run out of fuel and turn off. This allows 
the cooling flow to eventually resume and deposit a fresh supply of fuel to the AGN. 

In the absence of any disturbances (i.e., reheating of the X-ray gas), one would expect the 
gas deposited by a X-ray cooling flow to stream inwards along all radial directions. Instead, 
all three molecular filaments detected in Per A are confined to the east-west directions, 



reflecting the global distribution of the molecular gas in this galaxy (jSalome et al.ll2006l ). 
This is just the configuration expected if the cooling flow is quenched (or reduced below 
the level of detectability) over a large solid angle in the north-south direction by the radio 
jets from Per A, but is still able to operate in directions orthogonal to these jets. Indeed, as 
pointed out in §3.4, the bulk of the molecular gas detected coincides with a band of relatively 
cool X-ray gas lying between the X-ray cavities inflated by the radio jets. Our observations 
implicate this relatively cool X-ray gas as part of the cooling flow. 

Note that, when deprojected in three dimensions, the molecular gas is not infalling along 
all directions orthogonal to the radio jets, but is instead confined to certain radial directions 
mostly toward the far side of the galaxy. This suggests that, even orthogonal to the radio jets, 
large portions of the X-ray gas may be experiencing reheating, reducing their cooling flow if 
any to below the level of detectability. Indeed, the X-ray image of the Perseus cluster shows 
ripples in virtually all directions projected in the plane of the sky, but which is strongest 
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over a broad range of angles centered along the north-south directions (IFabian et al.ll2006l ). 
Each ripple is believed to be an acoustic wave driven by episodic expansions in the X-ray 
bubbles inflated by the radio jets from Per A. Once a given bubble reaches a sufficiently high 
pressure to become buoyant it rises through the X-ray gas, whereby the radio jet creates new 
X-ray bubbles at the center. The observed ripples (acoustic waves) are believed to gently 
heat the surrounding gas. Our results indicate that a (detectable) cooling flow is still able 
to operate in a few localized regions where presumably the X-ray gas is not being strongly 
reheated. Even at these locations, the cooling flow is likely intermittent or modulated in 
time by disturbances propagating through the X-ray gas. 



Dunn Sz Fabian! (120061 ) infer an average duty cycle for the bubbling activity in putative 
cooling- flow clusters of at least ~70%, and perhaps as high as ~90%. If the cooling flow is 
severely disrupted if not completely quenched when the AGN is on, such a high duty cycle 
may not permit the X-ray gas to both cool and reach the central supermassive black hole 
in the cD galaxy during the relatively short time when the AGN is off. Instead, our results 
imply that a residual cooling flow occurs even when the AGN is on, and that this cooling 
flow is able to provide a nearly continuously supply of fuel — most likely in the form of 
molecular gas — to the AGN. 



7. Summary and Conclusions 

To better understand the origin of the molecular gas in Per A, we have imaged this 
galaxy in CO(2-l) at an angular resolution of 3" (spatial resolution of 1 kpc) with the 
SubMillimeter Array (SMA). Single-dish maps with angular resolutions as high as 12" (4 kpc) 
have previously shown that the molecular gas is concentrated in the east-west direction, and 
extends ~40" (~14 kpc) in both directions from center. This gas shows no discernible ordered 
spatial-kinematic pattern, exhibiting blueshifted velocities both east and west of center, and 
redshifted velocities at or close to the center. Our SMA map (Figs. 1 and 2), which covers 
a region of radius ~28" (~10 kpc) centered on Per A, reveals that: 

1. the molecular gas extends to the edge of our field, and is concentrated primarily in 
three radial filaments with projected lengths ranging from at least 1.1-2.4 kpc spanning 
the galaxy center to 7.7 kpc. These filaments, along with several blobs, have a total 
mass in molecular gas of (4.2±0.2) x 10 9 M , with the three filaments alone containing 
~80% of the detected gas mass. 

2. the bulk of the molecular gas detected (i.e., all three filaments and all except two blobs 
of molecular gas) is aligned roughly east-west and lie between two X-ray cavities to the 
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north and south of center. This gas is aligned approximately orthogonal to the axis of 
the radio jets from the AGN in Per A. 

3. the two outer (western and eastern) filaments, and the other blobs of molecular gas, 
exhibit only blueshifted velocities. The velocities of both these filaments change linearly 
with radius to larger blueshifted velocities at smaller radii. Neither these filaments nor 
the blobs are coincident with detectable silhouette dust, suggesting that they lie at a 
relatively large inclination to the plane of the sky on the far (back) side of Per A. 

4. the inner filament exhibits both blueshifted and redshifted velocities that straddle the 
systemic velocity of the galaxy. Unlike the two outer filaments, the velocity of the 
inner filament does not change linearly with radius; furthermore, this filament exhibits 
smaller blueshifted or larger redshifted velocities at smaller radii. It is not clear whether 
this filament is coincident with silhouette dust associated with Per A. 

5. all the molecular gas detected in the east- west direction is spatially coincident with Ha 
gas, often lying against or close to local Ha brightenings. Where measured, specifically 
at the locations of the inner and western filaments, the velocities of the Ha gas are 
comparable with that of the molecular gas. The molecular gas also coincides spatially 
with a band of relatively cool X-ray gas lying between the two X-ray cavities. 

We interpret the results in the following manner: 

1. the filaments are unlikely to have been captured from one or more gas- rich galaxies. 
Gas captured in this manner should orbit and hence show a global rotation about the 
cD galaxy, as has hitherto been the case in all other elliptical galaxies that have been 
well mapped in CO. 

2. lying on the far (back) side of Per A, the two outer filaments are therefore flowing 
radially inwards towards the center of the galaxy. The only known process that can 
naturally deposit gas in this manner is a X-ray cooling flow. 

3. the velocity pattern of the two outer filaments can be easily reproduced as free fall in 
the gravitational potential of Per A, with their outer tips lying ~1 kpc downstream 
from where the cooling gas first decouples from the X-ray gas. Both filaments have a 
dynamical age of ~20 Myrs, implying a total mass-deposition rate into these filaments 
of roughly ~75 M Q yr -1 (likely accurate to only a factor of a few). 

4. the velocity pattern of the inner filament cannot be explained with a simple model. We 
speculate that this filament likely traces gas settling into the gravitational potential 
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well of Per A, and perhaps corresponds (in part) to the outer regions of a putative hot 
molecular hydrogen gas disk with a radius of ~50 pc around the AGN. 

Despite the strong and widespread heating of the surrounding gas inferred from X-ray 
observations, our results provide the most direct evidence yet for a (residual) X-ray cooling 
flow in the Perseus cluster. This cooling flow may provide a nearly continuous supply of fuel 
for the AGN in Per A, which in turn regulates the cooling flow. 
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image, and C. J. Conselice for a fits file of the Ha image, used in this paper. We thank the 
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a visit to the Harvard CfA, and acknowledges a grant from the National Science Council of 
Taiwan in support of this work. The Submillimeter Array (SMA) is a collaborative project 
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Table 1. Parameters of Molecular Features 



Name 


Inner radius 

(asec) 


Outer radius 
(asec) 


Position angle 
(degrees) 


Inner velocity 
(km s- 1 ) 


Outer velocity 
(km s- 1 ) 


Flux density 
(Jy km s -1 ) 


Gass mass 
(xlO 8 M Q ) 


Inner filament 


1.2 ±0.3 


4.2 ±0.3 


~270 


+100 


-100 


62.1 ± 3.6 


18.0 ± 1.1 


Western filament 


9.0 ±0.3 


15.7 ±0.4 


~281 


-180 


-80 


29.3 ±2.6 


8.5 ± 0.8 


Wl 


21.6 ±0.6 




~198 


-20 




9.9 ± 1.1 


2.9 ± 0.3 


Eastern filament 


15.2 ±0.4 


21.3 ±0.3 


~76 


-160 


-40 


21.0 ±2.6 


6.1 ± 0.8 


El 


25.7 ±0.5 




~76 


-70 




14.9 ± 1.7 


4.3 ± 0.5 


E2 


22.8 ±0.7 




~63 


-80 




4.2 ±0.8 


1.2 ± 0.2 


Nl 


18.2 ±0.4 




~358 


-30 




4.9 ±0.8 


1.4 ±0.2 


SI 


20.9 ±0.6 




~211 


-130 




8.4 ±1.3 


2.4 ±0.4 



Note. — The flux densities listed have been corrected for the primary beam response of the SMA. The inner/outer radii and velocities 
listed for the filaments are derived from the channel maps. For the features Wl, El, E2, Nl, and SI, which are mostly spatially unresolved, 
we list only the radii and velocity at their centroids as derived from the moment maps. 
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Fig. 1. — Channel maps showing contours of CO(2-l) intensity with the velocity in- 
dicated in e ach panel measured with respe ct to the systemic heliocentric velocity of 
5264 km s -1 (IHuchra. Vogeley. fc Gellerl 119991 ). Contour levels are plotted at -3, 3, 4, 5, 
6, 7, 8 x 16 mJy beam -1 (root-mean-square noise level; la). No primary beam correction 
has been applied to these maps. The gray background corresponds to the integrated CO(2-l) 
intensity after primary beam correction as shown in Figure 2. The position of the nuclear 
continuum source in Per A is indicated by a cross. The synthesized beam is shown as an 
ellipse at the lower left corner of the top left panel, and has a size of 3'.'0 x 2'!7 and position 
angle 39°. 5. The filaments indicated in Figure 2 are not uniformly bright in intensity along 
their lengths, and hence may not appear continuous in these channel maps. 
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Fig. 2. — Contours of the integrated C0(2-l) intensity plotted at levels of 3, 6, 9, 12, 15, 
20, 30, 40, and 50 x 0.32 Jy km s -1 , and color-coded map of the intensity-weighted C0(2-l) 
mean velocity measured with resp ect to the systemic heliocentric velocity of 5264±11 km s _1 
(jHuchra. Vogeley. fc Geilerlll999l ) as indicated by the right vertical bar. The full-width half- 
maximum of the SMA primary beam is 55" (20 kpc), as indicated by the dotted circle. The 
integrated CO (2-1) intensity has been corrected for the primary beam response. The map is 
centered on the active nucleus of Per A (Right Ascension = 03 h 19 m 48 s . 16, Declination = 
41° 30' 42'.'1), which was detected in the continuum at the location indicated by a cross. The 
synthesized beam is shown as a filled ellipse at the lower left corner, and has size of 3'.'0 x 2'/7 
with major axis at a position angle of 39°. 5. The inner, western, and eastern filaments, along 
with the features El, E2, Wl, Nl, and SI, are labeled as described in the text. 
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Fig. 3. — Spectra of all the features labeled in Figure 2, derived from the channel maps 
shown in Figure 1 but after correcting for the primary beam response of the SMA. An error 
bar corresponding to ±ler uncertainty is plotted in each panel. 
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Fig. 4. — Position-velocity (PV) diagrams showing CO(2-l) intensity along the radial extent 
of (a) inner filament along a position angle of 270°, (b) western filament along a position 
angle of 281°, and (c) eastern filament along a position angle of 76° as a function of velocity. 
Contour levels are plotted at —2, 2, 3, 4, 5, 6, and 7 x 16 mJy beam -1 (la). Positions 
correspond to offsets from the center of Per A with negative values indicating west, and 
yelocities are measured with resp ect to the systemic heliocentric velocity of 5264±11 km s _1 
( iHuchra. Vogeley. fc Geilerlll999l ). In panel (b), the feature at positions > —5" and velocities 
> —100 km s~ x corresponds to the inner filament. In panel (c), the feature at a position 
of 26" and velocity of —80 km s _1 that is discontinous in velocity from the eastern filament 
corresponds to El (Fig. 2), which is probably a separate structure as described in the text. 
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Fig. 5. — Contours of the integrated CO(2-l) intensity (from Fig. 2) overlaid o n a fal se-color 
map of the Ha nebula associated with Per A adapted from IConselice et al.l (120011 ). This 
maps spans only the inner central region of the Ha nebula. 
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Fig. 6. — Black contours of the integrated C0(2-l) intensity ( from Fig. 2) and red contours 
of the continuum emission at 90 cm from iPedlar et al.l (119901 ) s howing radio jets fro m the 
AGN in Per A, both overlaid on a false-color X-ray image from iFabian et al.l (120061 ). The 
synthesized beams of the CO (2-1) and 90 cm images are shown as the black and red ellipses, 
respectively, at the lower left corner. 
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Fig. 7. — Contours of the integrated CO(2-l) intensity (from Fig. 2) overlaid on a X-ray 
map of the inferred mass distribution of the 0.5 keV (~5 x 10 6 K) component in the core of 
the Perseus cluster (from lFabian et al.ll2006l ). 



Fig. 8. — Contours of the integrated CO(2-l) intens i ty (fr om Fig. 2) overlaid on the HST- 
WFPC2 optical image of Per A from ICarlson et al.l (119981 ). Dark features are dust seen in 
silhouette against the optical continuum (primarily starlight) from Per A, with a possible 
contribution also to the optical continuum from a forego und infalling galaxy (the HVS; see 
text). The blob SI lies outside the region spanned by this optical image. 
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Fig. 9. — Contours of the integrated CO(2-l) intensity (from Fig. 2) overlaid on the 
HST-WFPC2 B-ban d transmission map (see text) of the central region of Per A from 
Keel &: White! ( 120011 ) . Bright features are intensities above a smooth background, and dark 
features below the same background. Only the inner and western molecular filaments, along 
with the blobs Wl and Nl, are visible over the region spanned by this optical map. 
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Fig. 10. — Velocity of an object in free-fall along the gravitational potential of Per A as 
a function of radial distance from center, computed from Eq. (4) for two different galaxy 
masses as indicated by the solid and dashed curves. Panel (a) shows the computed velocity 
profile of an object dropped from a radius of 8.5 kpc, and superposed the PV-diagram of 
the western filament (contours; from Fig. 4) deprojected by an inclination of ~40° so that it 
spans 4.2-7.4 kpc as indicated by the double-headed arrow. Panel (b) shows the computed 
velocity profile of an object dropped from a radius of 12.2 kpc, and superposed the PV- 
diagram (contours; from Fig. 4) of the eastern filament deprojected by an inclination of 
~47° so that it spans 8.1-11.3 kpc as indicated by the double-headed arrow. As can be 
seen, the computed velocity profile assuming a lower galaxy mass of 3.4 x 10 11 M Q (solid 
curve) better matches the observed PV-diagrams than the computed velocity profile for a 
galaxy mass of 8.3 x 10 11 M inferred from the luminosity of the galaxy and its estimated 
mass-to-light ratio as described in the text. 



